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1. Introduction. 

The main difficulty connected with the attempt to fly by means of a 
machine heavier than air is that of longitudinal stability. It is not 
difficult to construct an aeroplane system which shall be transversely 
stable. 

But for this difficulty the problem of artificial flight would probably 
have been solved already. Experiments in gliding under gravity have 
been always made with machines not too large to be kept balanced 
by the skill of the experimenter, and the glides, though undoubtedly 
successful, have been, of short duration. Experiments have invariably 
stopped short of the performance of continuous flight by a mechanically 
propelled machine. 

The problem of artificial flight is hardly likely to be solved until the 
conditions of longitudinal stability of an aeroplane system have been 
reduced to a matter of pure mathematical calculation. 

A theoretical investigation, even if calculated under conditions 
slightly different to those occurring in nature, will serve as a basis 
of comparison by which experimental results can be co-ordinated and 
interpreted in their true light. 

The object of these investigations is (1) to show that the longitudinal 
stability of aeroplane systems can be made the subject of mathematical 
calculation; (2) to draw the attention of those interested in the 
problem of artificial flight to the necessity of acquiring further 
experimental knowledge concerning the quantities on which this 
stability is shown to depend. 

2. General Investigation of the Longitudinal Stability of any Symmetrical 

A erro])lane System, 

Consider any system of aeroplanes, having a plane of symmetry, 
descending in a vertical plane, in air or in any resisting medium 
whatever. To specify the motion (which we suppose to be two- 
dimensional) choose two axes at right angles fixed in the body, having 
the centre of gravity as origin. 

The motion will be completely determined if at every instant we 
know— 

(1) the angle which the axis of x makes with the vertical.* 

(2) the velocity components %, v, of the body along the two axes. 

* This angle is supposed measured from the dotvnwarcl drawn vertical in the 
positive direction. The axis of y must be drawn above the horizontal at an 
inclination of #. 
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The angular velocity of rotation of the body will be dO/dt or 6. 

We shall use m to denote the mass of the body, mk 2 its moment of 
inertia about the centre of gravity. 

Whatever be the law of resistance, the resistances of the air on the 
several parts of the system will in general be functions of u, v and 6. 
These resistances are always equivalent to two forces, which we shall 
call mX and mY, along the axes, and a couple mG about the origin, so 
that X, Y, G, denote the forces and couple, divided by the mass of 
the body. 

The equations of motion of the body are — 

(du d0\ n v 

m -77 - v -rr) = mg cos - mX, 
\dt dtj y 



!dv dO 
m 



+ n-r, ) = - niq sin - mY. 
dtj J ' 



7 o U/ U ft 

mk 1 -— = - mG. 



> •*••*•»»« f X /• 



In steady motion u, v and are constant, and equations (1) give 

= #cos<9-X, = ~-#sin<9~Y, = -G (2). 

Knowing the forms of the aeroplanes and other parts of the system 
and the law of resistance, X, Y, G are known functions of u 9 v and 0, 
Moreover, in steady motion, = 0. 

Equations (2) thus determine the values of u, v, 0, for steady 
motion. 

Fluctuations about Steady Motion. — We must now examine what 
happens when the system is slightly disturbed from its state of steady 
motion. Let the disturbance be represented at time t by small 
increases 8u, &», 86 in the values of u, v, (9, and let %, v , be their 
values in steady motion, so that in the disturbed motion, 

u — u + 8u, v = v + 8v, 6 = + 80. 

Also let X , Y , G be the forces and couple in the steady state ; 
then in the disturbed state we have, neglecting small quantities of the 
second order, 

X = X + &u— + 8v ( ^ + 8e (j ^ (3), 

du do dO 

and two similar equations for Y and G. 

We shall denote differential coefficients such as dX./du . . . dQ/d$ 
by jx-u . . . lx# . 

Substituting in the equations of motion, we obtain, to the first 
order, 
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--- 8u - v -1 80 = - sin # S# - X u Sw - X v 8v - Xe80, 
at at 

§L 8v + uo i-80 = -g cos 6 86 - Y^ - Y v 8v - Ye$0, 
at at 



» I <"> 



G %l 8u - G v 8v - Gre86. 



To solve these equations, put 

8u = Ye* , &> = Q € *< , 8(9 - R € *' . 

Substituting, dividing by e xt , and re-arranging the terms, we have 

P (X + X u ) + QX-y + R (XX# - X% + # sin O ) = 0, 

PY tt + Q (X + Y„) + R (XYe + X^ + gr cos <9 ) = 0, 

PGu+QG^ + K^ + XG*) = 0. 

Eliminating P, Q, R, 

X + X tt , X.y , - Xv + XXg + g sin # 
Y M , X + Y v , X2^ + XYg + # cos # 
G^ ? Gy , X 2 & 2 + XG# 



= (4). 



If the determinant be expanded in powers of X we get an equation 

of the form 

AX4 + BX3 + CX 3 + DX + E - 



«•«*«« 



o «* » s » 



(4«), 



where 



A = F, 

B = k*(Xu + Y v ) + Ge, 

C = ¥ (X U Y V - X,Y„,) + * G„ - %G, - X* G tt - Y^ G„ 

D = u Q (X v G u - X U G V ) + v (Y v & u . - Y„G,>), J. (5). 

- g sin G U - g cos (9 G V + Ye (X V G U - X«Gy) 

— X# (Y V G U — Y U G V ) + Gg (X u Yv - X V Y W ), 
E = £ cos #o (X v G tt - X tt G,>) - g sin <9 (Y V G U - Y„G„), j 

In many cases it is possible to take the axes of co-ordinates, so that 
shall be zero, and also that Y and its differentials shall vanish. In 
such cases the coefficients take the following simple forms : — 



A = k\ 

B = WX U + Ge , 

C = v G u - uqG v - X# G m { 

I) = u (X V G U - X U G ) - gG m 

E = o {XuG u — j\ u G u ), 



6). 



1904.] The Longitudinal Stability of Aerial Gliders. 103 

Equation (4) or (4a) is the period equation for small fluctuations 
about steady motion. 

In order that the steady motion may be stable, the roots of this 
biquadratic must either be real and negative, or complex with their 
real parts negative, and Bouth* finds that this will be the case if the 
six quantities 

A, B, C, D, E, and BCD - AD 2 - EB*, 

are all of the same sign. Since A is essentially positive the remaining 
five quantities must all be positive. We shall denote the last quantity 
byH. 

3. General Theorems A 

(1) The general transformation formulae, connecting the nine co- 
efficients X tt . . . G<j referred to any given system of axes with those 
referred to any other system of axes, may be easily written down and 
need not be discussed here. 

(2) The work of calculating these coefficients for a given system 
may be reduced if the resistance is proportional to the square of the 
velocity, for X, Y, G- will then be homogeneous quadratic functions of 
u, v, v, and Euler's theorem of homogeneous functions gives, remember- 
ing that $o = and applying (2), 

iiX u + vX v = 2g cos By 

uY u -hvY v =- -2g sin #, 

uG M + vQ v = 0. 

(3) If V is the velocity of gliding, the coefficients mK u ... mGe are 
all linear functions of V for a given angle of gliding. But mg being, 
in steady motion, equal to the vertical resistance, is proportional to V 2 . 
Hence if m be eliminated, the values of X M . . . G<? are inversely 
proportional to V. In this case 

A is independent of V, 

is of dimensions V"" 1 , 

is of form Jr + Cjv , 
D is of dimensions V™ 1 , 
E is of dimensions V" 2 , 

and the expression H or BCD - AD 2 - EB 2 assumes the form 
py-2 _f_ QV~ 4 . Thus the conditions of stability C>OandH>0 
impose limits on the value of V 2 , the remaining conditions only 

* Routh, ' Advanced Rigid Dynamics,' p. 167. 

f This section and also the subsequent parts enclosed in [.;....] have been 
rewritten October 27, 1903.° 
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depend on the form and dimensions of the machine and the angle of 
gliding. 

We shall now show how the nine coefficients X w . . . G# can be 
calculated for a system of aeroplanes if the laws of variation of the 
resultant pressures and of the positions of the centres of pressure are 
known. We shall assume that the resultant pressure on a plane 
lamina varies as the square of the velocity, and acts in a direction 
normal to the lamina. 

We may therefore write 

E - KSV 2 /(a), 

where S is the area of the lamina, V the velocity, and a the angle 
between the plane of the lamina and the direction of motion, K a 
constant depending on the units employed. 

The function / (a) has been determined by Langley* for certain 
rectangular planes. 

Also, let the distance of the centre of pressure from the centre 
of figure be a<fi (a), 2a being the breadth of the lamina. 

Experiments to determine <f>. ( a ) fo f square planes have been made 
by Joesel,t Kummer,{ and Langley,§ and Kummer has also experi- 
mented on oblong planes, with their longer side in the direction of 
motion. 

Their results show a certain amount of discrepancy, and there 
appears to be considerable difficulty in obtaining consistent results 
at small inclinations. 

Further experiments are very necessary, and it is to be hoped that 
more attention will be given to determinations of <j> (a), when their 
importance, as affecting the stability, has been recognised. 

Experiments show that / (a), <£ (a) are, to a first approximation, 
independent of the translational velocity of the lamina, but the effects 
of a rotational angular velocity, 6, have never been considered. 
It would not be difficult to determine these effects by experiments 
with a whirling table, making 6 the angular velocity of the table. 
Eailing such experiments, these effects must be neglected, and, as 6 is 
jzero in the steady motion, and only small oscillations are considered, 
they are probably small. When the system consists of a number of 
planes rigidly connected together, we may now write the component 
forces and couple in the form 

mX = 2KSiVi 2 /(ai - ft) cos ft ^ 

mY = SKSi Vi*/ (ai - ft) sin ft I ..... , (7), 

mQ = 2KSiVi 2 /(ai - ft) {p x + a x <f> (a x - ft)} ^ 

* ' Experiments in Aerodynamics,' p. 62, 
f Joesel, ' Memorial du Genie Maritime,' 1870. 
% * Berlin Akad. Abhandlungen/ 1875, 1876. 
§ Zoc. cii,, p. 90. 
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where Si is the area, 2&i the breadth of the lamina, £1,1)1 the co- 
ordinates of its centre, Vi the resultant velocity of its centre, 
a\ — tan™ 1 Ui/v-i the angle which the direction of motion makes with 
the axis of y, ft the angle which the plane of the lamina makes with 
the axis of y, so that (&i - fa) is the angle between the plane and 
the direction of motion of the centre of mass, and finally, 
2h = 171 cos fa + £1 sin fa. 

The summation is to be extended over all the planes of the 
system. 

If, now, u, v be the component velocities of the centre of gravity 
the velocities it u Vi of the centre of the lamina will be given by 

and for steady motion, since = 0, u\ = u, V\ — v. 

We may therefore write the expressions for the forces in the form 

mX = 2K$i (tti 2 + v 1 2 )f(t f <m~~ 1 u 1 /v 1 - fa) cos fa, 

= 2KSi (^ + ^-2^+2^// tan-if!^Il3^ -ftlcosft, 

L ^ + ^/ J 

neglecting terms in # 2 . 

Similarly, 

•wY = 2KS X (# + # _ 2w7?1 d + 2^id)/ ( tan- 1 l^LlJ} ) -ft "I sin ft , 

L W + ^ft/ J 

* 

<mG - - 2KS 2 (# + & - 2u Vl + 2u^0)f ( tan- 1 (^J.f) - ft ) 

L \ v+ t l oJ J 



{»♦*(<-■* &?!)-*■)}• 



■ # + £i#' 

Differentiating these expressions with respect to %, v, we have :- — 
mX u = 2 [KSi (2« - r h 0)f(a 1 - fa) cos ft + KSi%/' (a 2 - ft) cos ft] ; 

•or, since = in the steady motion, 

mX^ = 2 [2KSi%/(ai - ft) cos ft + KS#/ , (a 1 - ft) cos ft].... (S), 

Similarly, 

■mX v - 2 [2KS i «;/(a 1 - ft) cos ft - KSi«/'(ai - ft) cos ft], 

mY tt = 2 [2KSiw/(ai - ft) sin ft + KSi^/"' (a x - ft) sin ft], 

.mYtf = 2 [2KSiv/(a! - ft) sin ft - KBiuf (gc 1 - ft) sin ft], 

mXe = 2 [2KSi ( - wgi + v^ 1 )f(a 1 - ft) cos ft 

- KS X («?J7i + w£i)/'(«i - ft) cos ft], 
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mY- e = 2 [2KSi ( - urn + v^)ffa - ft) sin ft 

- KSi (^i + ^i) /' fa ™ ft) sin ft], 
mQ u = 2 [ - 2KSi (% - yj)ffa - ft){>i + %</> (ai - ft)} 

- KSi^/' (ai - ft){f i + %<£ (ai - ft)} 

- KSviffa - ft) a^'fa ™ A)]- 

Since in the steady motion G- = 0, 

£<?., 2 (u 2 + 'y 2 ) 2KSi/(a! - ft) {pi + a-i<t> fa - ft)} = 0, 

mQc u = #2 [ - KSi/' (a! - ft) {^ + a^ (ai - ft)} 

- KSi/(ai - ft)ai<^'(ai - ft], 

wiG* - %S [KSi/' (ai - ft.) {jpi + %<£ (a! - ft)} 

+ KS 1 /(« 1 -ft)a 1 *> 1 -ft)], 

wiG* = 2 [ - KSi /' (a! - ft) ( - ?% + «£i) {p + acf> fa - ft)} 

- KSi/(a! - ft) ( - ^x + %&) aft fa - ft)]. 

4. Numerical Calculations for Particular Cases, Single Lamince. 

We now consider the stability of certain particular systems, 
beginning with a single plane lamina. 

Since the lamina is falling steadily under gravity, its plar:3 must 
necessarily be horizontal. 

Take the axis of y parallel to this plane, then = 0, /3 = 0, and, 
therefore, Y — ; also, for equilibrium, rj + a<j) (a) = 0. 

The coefficients, therefore, become 

mXu = 2KBuf(a) + KSvf(a) i 

mX v = 2KSvf fa) - KSuf (a), 

mXg « - 2KS (w»j - v£)f(a) + KS/'(a) (wy + t*£), 

mG u = - KS#/(a) a <f> (a), 

mGy = KBuffaafifai 

mQc'e = -KSa/(a)(w7 + w£)a<£'(a). 

iste. 1.— Consider a square plane, balanced so as to fly at an angle 
of 10°, with the centre of gravity in the plane, and suppose its radius 
of gyration given by k 2 = ^a 2 . 

Professor Langiey's results give for this particular angle 

/(a) = 0-3, /'(a) = 1-6, 

and from JoesePs formula 

<£(*) - 0-49, <£'(a) = -0*59, .\ 17 = 0'49a. 
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also 

u = V sin a = 0-17V, v = 0'98V. 

Substituting these values in the above expressions for X„ . . ., we 

mX u = 1-66KSV, mG tt = 0'l7KSYa, 

mX„ = 0-32KSV, mG„ - -0-03KSYa, 

mX* - 0*73KSYa, mG* = 0«09KSVa 2 . 

Substituting these values in the expressions for A, B, C, D, E, given 
in § (2), and remembering that by the conditions for steady motion 
mg = KSY 2 /(a), we obtain 

A = Ja 2 , B = 92 a 2 /Y, C = 18a - 1340 a 2 /Y 2 , 

D = 295 a/Y, E = 36400 a/Y 2 , 

and therefore 

H = BCD - AD 2 - EB 2 = 4-8 . 10 5 ~ - 3-6 . 10 7 ^ 

-4-3. 10 4 iL- 3-08. 10 8 ~. 

This expression will be positive if V 2 > 774a, and this condition will 
also make C positive. The glider will therefore be stable if its 
velocity is greater than ^(7 74a), the units being feet and seconds. 

Ex. 2. — Let us now take the angle of gliding to be 35°, and assume as 
before k 2 — -| a 2 . At this angle we have 

/(a) = 0-84, /'(a) =0-6, 

<£ (a) - 0'26 ? </>' (a) = 0'5. 

With these values, we have 

mX u = 1-34KSV, mG u = 0'35KSYa, 

mX, = MKSY, mG, - 0-21KSVa, 

mXe - -0-lKSYa, mG* = 0'08KSYa. 

Substituting in the expressions for A, B, etc., we have 

A = |-a 2 , B = 27 a 2 /V, 

C = 15-2a + 50 a 2 / V 2 , D = 380 a/Y + 240 a/Y = 620 a/Y, 

E = 25600 a 2 /Y 2 , H - 5600 a 3 /Y 2 -- 17400000 a 5 /Y 4 . 

The latter is positive when Y 2 > 3100a, which is therefore the con- 
dition of stability. 

Ex. 3. — Let us now take an oblong plane lamina ; the values of / (a) 
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are given by Langley for a lamina of 30" by 4/8" the shorter side 
being in the plane of motion. In the absence of any information as to 
<£ (oc), we shall assume (whether correctly or incorrectly) the same 
values as for square planes. 

Let the moment of inertia be given, as before, by k 2 = \a 2 . 

For an angle of 10° we have 

f(oc) = 0-44, /» = 1-5, 

<£(a) = 0-49, <£'(a) = 0'59. 

Substituting these values in the expressions for the coefficients, we 
have 

mX u = 1-6KSV, mG u = 0-26KSVa, 

mX v = 0-6KSV, mG v = - 0-4KSVa, 

mXe = 0'56KSVa, mG« = 0'13KSVa 2 . 

Substituting again in the expressions for A, B, C, D, and E, we 
have, assuming k 2 == |a 2 , 

A - F = i« 2 , 

B - k 2 X u + G' e = I(l-6~ + 0-13a 2 \KSV = 60a 2 /V, 

C = flG tt -%G„-X* G w = 19 -4a -7 60a 2 / V 2 , 
D = u(X v G u -X u G v )-gG v = 280a/V, 
E = # (X V G U - X^G,) = 35000a/V 2 . 

Therefore 

H = 3-2. 10 5 ^-l-2.10 7 ^ 4 -3-9. 10 4 ^-1'2. 10 8 ^. 

This is positive if Y 2 > 470«, which is therefore the condition of 
stability. 

5. Gliders Formed of two Planes. 

Proceeding now to consider the stability of systems made up of 
several planes rigidly connected together, we shall first consider the 
stability of a gliding system supported on two slats Si S 2 , which are 
so narrow that displacements of their centre of pressures due to varia- 
tions of the angles of incidence of the wind may be neglected. 

^1 V2 



We shall first suppose that the two slats are in the same plane, and 
that the centre of gravity is also in this plane. 
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Putting <fi (a) = and ft = ft 2 = in the expressions of § (3), we 
have all the Y's equal to zero and 

mX u = K(Si + S 2 )(2^/(a) + fl/'(a)), (9) 

mX, - K (Si + S 2 ) (2*;/ (a) + < (a% 

mX§ = - K (Si^i + S 2 ^ 2 ) (2«/a 4- #/'a), 

*« = - W' (a) (S^i + S 2 i? 2 ), 

mG* = KtA/'(a)(Sir7i + S 2 ^2), 
mG* = Ka/' (cc) (S^i 2 + S 2 ^ 2 2 ). 

Now for equilibrium we must have 

Si^i-fS 2 ^ 2 = 0. 

In this case, G w , G^, and X§ vanish, and therefore, from (6), we see 
that the coefficients C, D, E, will be zero, and the equilibrum is critical 
or neutral. 

[This result, which is also evident from first principles, holds good 
equally when the centre of gravity is not in the plane of the laminge. 
In order to make the system really stable, the laminas must be inclined 
at small angles to the line joining their centres. 

Calling these angles ft and ft, neglecting <j> (a), and writing 
a' == a - ft, a ' = a - ft, we get in the equations (8), 

mX u = K2S X (2uf(cc f ) + vf(a!)) cos ft, 
mX v = K2Si (2vf(a!) - tif (a/)) cos ft, 
mX* = -K2Si^i(2i*/(a / )+^/'{a / ))cosft, 
mY tt = K2Si (2uf(af) + #/' (a')) sin ft, 
wY„ = K2Si(2^(a')-^/'(a'))sinft, 
mYi = - K2Si^x (2uf(af) + vf(a')) sin ft, 
mG tt = -KvSSipi/^a'), 

mGd =* K^SSi^ijpi/^a'). 
The conditions for steady motion give 

= kSrf (a) sin fo + kS 2 f (a") sin f3 2i 
mg[V* = ^SV(^)cosft + /^ 2 /(^)cosft, 

= &Si/(a') Tyx cos /3 + kSzffa") rj 2 cos ft. 

With these substitutions, those of the nine coefficients X tt . » . G# which 
vanish in the limiting case of coplanar laminae are given by 



x# 


Qu - 


V 


V 


-*~U 


„ Y * 


V 


u 
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G, = _ g jm^ //'(*') _/'(*") 

# sin ft sin ft f /'( a ') _/'( a ") 

and an expression for Y# , which is not required. 

The remaining coefficients may for a first approximation be taken 
the same as in (9). We notice the following points. 

(1) The expressions for X# f Y u , Y v> Gr u and Gt v contain the factor 

/'(a') //(a") — /'(a")//( a ")» ^ n the -^ aws °^ resistance commonly 
assumed /'(a) //(a) decreases as a increases, and this coefficient, 
therefore, becomes negative if a! < a'', i.g., if the two planes slope towards 
each other a! referring to the front plane. This is the case in a gliding 
machine furnished with a rudder, if the rudder is tilted slightly 
upwards. In such cases t/ 2 must be negative, hence X# and G u are 
positive, Y u and G u are negative. 

(2) If ft and ft are of the first order of small quantities, ft - ft 2 
also being of the first order, the expressions X$ f Q u and G v are of the 
first order, but Y u and Y^ are of the second order, and hence equations 
(6) give approximately the five coefficients A, B, C, D, E. 

(3) The fifth stability condition H > O appears at first sight difficult 
to reconcile with the smallness of the coefficients C, D, E for small 
inclinations of the planes, as B 2 E is of a lower order of small quantities 
than BCD. To satisfy this fifth condition, however, the important 
thing is to make E small compared with C. Now C will be found to 
consist of two parts, one independent of Y and the other negative and 
proportional to 1/Y 2 , while E is proportional to 1/Y 2 , hence stability 
can best be secured by making Y 2 sufficiently large.] 

(4) If the planes are not infinitely narrow, it will be found that 
X tt , X y , Y u , Y v are the same as before, and - mQ u /v, and mG v /u are 
both increased by 2KSi«i (/'( a ') <£ ( a ') + /(a') <£' (&))> This will have 
the effect of altering C, E, D in the same proportion, and will, there- 
fore, not alter the stability, except when a = 0, when the additional 
term will prevent E]D from becoming zero. Ge is increased by the 
term #2KSi<2i?7i (/'(#')<£(#') -t-/(a')<£'(a )). This term will, in general, 
be negative, and therefore GU will be diminished, which will have the 
effect of diminishing B. 

6. Examples of Two-Plane Gliders. 

We shall now assume certain particular values for the dimensions of 
gliders of this form, and proceed to calculate the conditions of 
stability. 

Ex. 4. — Let us first consider the case of the two slats of equal area 
set at angles such that 

II-' O f- o 

5, cc 2 = 5 . 
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Then we have from Langley*s results, 

/(*') = 0-56, /'(«') = 1, 

f(a! r ) = 0*28, /'(a") = 2-5; 

for equilibrium 

a'/ (a") - - *?2/0?i) ; 

therefore, if we put ?}i = a, we have ^ = 2a. Let & 2 be assumed = 2a 2 . 
Substituting in the expressions for X w etc, we have 

mX„ = 3*6KSY, mG w - 4KSYa, 

mX, = 1KSY, mG, = 0-7KSYa, 

wXfl » 4KSYa, roG* = llKSYa 2 , 
and therefore, 

A = 2a 2 , B = 680 ^, C = 156a - 22800a 2 / V 2 , 

D-2721* E = 297000 ~ 9 , H = 2,7 . 10»£ >- 1-4. 10* .l* 

V V 2 V 2 Y 2 

-2'9.10 10 ^-l-3;10 u £, 

This is positive if Y 2 > 590a. 

Ex. 5. — Let us now take the case of two equal square planes inclined 
at a small angle to each other. 

Let a = 10°, ft = -5°, ft = 5°, 

and let 2a be the breadth of either plane. 
For these angles we have 

f(«i~Pi) = 0-U, . />i~-ft) = 1*7, 
/(«2-)82) = 0-15, /'(«2-ft)-l'8, 

^ ( ai - ft) = 0*45, £' ( ai - ft) « 0*6, 

<^> (as - ft) - 0-55, <£' ( ai - ft) = 0-58. 

For equilibrium we must have 

/(ai - ft) {171.+ a<£ (ai - ft)} = f(a 2 - ft) {^ 2 + a</> (oc 2 - ft)} , 

whence 

3^1 + ^2 = 0'28a; 

therefore, if we put v% = - 4a, we have ^i =. 1 *4a. Let & 2 be assumed 
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to be = 4a 2 . Substituting these values in the expressions for the 
coefficients, we have 

mX u - 3-6KSV, mY u = 0-02KSV, 

mX v = 0-58KSV, mY v = -O'lKSV, 

mXe - l'SKSVa, mYe - 0-OOKSV, 

mGt u = 2-61KSVa, mG, = 0'44KSVa, 
mG* = 30KSVa 2 . 

Substituting again in the expressions for A, B, C, D, E, we obtain, 
putting h 2 = 4a 2 , 

A - 4a 2 , B = 2390 ~, C = 142a- 13700 |i, 

I) - 2500^, E = 284000^,. 

H is positive if V 2 > 2000a, which is the condition of stability for 
this form of glider. 

By putting I = 7*4a, Z represents approximately the extreme length 
of the glider, and the condition of stability reduces to V 2 > 2701, 

Ex. 6. — Consider next the case of two unequal square planes, 
inclined at an angle, and in the first case suppose that the smaller 
plane is in front. 

Let S and S' be the respective areas, and suppose S = 10S f . 

Let 2a be the breadth of the large plane, and 2a/ 3*1 that of the 
smaller, also let the distance between the centres be equal to 3a : and 
let Jc 2 — a 2 . 

Let the angle a between the large plane and the direction of motion 
be 10°, and let the angle between the two planes be also 10°, so that 
the small plane is inclined at an angle of 20° to the direction of 
motion. 

Then we have 

f(a) = 0*3, /'(a) = 1-6, <£ (a) = 0*49, <£'(a) = 0*59, 

/(a + 0) = O-5, /'(* + /?) = 1-57, <£(a + /3) = 0'4, <£'(a + /3) = 0*56. 

For equilibrium we have, if ?? be the distance between e.g. and centre 
of large plane, 

S/(a) ( v ~ a$(oc)) - S'/ (a + (*)(b + a' <j>(a + /3)- v ). 

Therefore >? = # 94a, so that ^ = i - rj = 2*06a. 
Substituting in the expressions for X w , etc., 

mX n = 1-87 KSV, mX v = 0-41 KSV, mX* = 1-4 KSVa, 
mG„ = 0-57 KSVa, mG„ = 0-097 KSVa, mG* = 1*13 KSVa 2 . 
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Substituting these values in the expressions for A, B, etc, assuming 
that k 2 = a 2 we have 



ur „ _ a? 



A = a 2 , B = 273~, C = 50-8a- 15050 ys 



? 



D = 917--, E = 120000^, H - 118.10^-123.10^. 

This will be positive if V 2 > 1040a. This condition will also make 
C positive, and is therefore the condition of stability. 

Ex. 7. —Let us now suppose the Small plane to be placed behind the 
other 

s $' 



being inclined to it at an angle of 5°, the direction of motion again 
making an angle of 10° with the large plane and k 2 being = a 2 , as 
in Ex. 6. 

In this case we have 

/(a) = 0-3, /'(a) = 1*6, <£(a) - 0*49, <//(a) = 0"59, 
f(a') - 0-17, /'(a') = 1-8, </>«> = 0*55, (//(a') - 0'6. 

If the distance between the centres of the planes be* 3a, the condi- 
tions of equilibrium give 

t) = -0*35a, rjj = ~2*65a. 

Substituting, we obtain 

ywX w = 1 -9 KSVa, mG tt - 0-33 KSVa, 

wX, - 0-41 KSVa, mG i} = 0"056 KSVa, 

wX -0-65 KSVa, mG„ - 1'32 KSVa 2 . 

,\ A = a 2 , B^233^, C = 34'8a - 2400^, 

V V 

D = 1780 " E = 80000 " , H = 12 . 10« f - 3 . 10 9 ^ . 

V V^ V 2 V 4 

H will be positive if V 2 > 250a. 

The condition of stability is therefore that V > x /(250a). 

[7. Effect of Moment of Inertia on Stability. 

It will be seen from (5) that the radius of gyration, k 2 , occurs only 
in the expressions for A, B, the first two coefficients of the deter- 
minantal equation. We have A = / t ; 2 , B = k 2 (X u + Y v ) + Ge. In all 
the cases considered X„ + Y v and G d ' are positive, and therefore these 
two coefficients are positive for all values of fc 2 . 
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Taking the expression 

H - BCD - AD 2 - EB 2 , 
it may be written 

CDGe - EG* 3 + *2 i(X tt + Y c ) CD - D* - 2 (X H + Y„)G*E} - **(Xu + Y^E. 

H is therefore increased with increase of ¥ if 

(X* + Y v ) (CD - 2GiE) - D 2 - 2fc 2 (X w + Y v ) E > 0. 

At the critical velocity, if II = 0, this becomes 

GDQb + (X„ + Y r ) 2 E/^ < EG* 2 . 

Unless this condition is fulfilled, the critical velocity given by H - 
increases with increase of k 2 . 

In all the numerical examples considered above, CDG* >EGr and 
E is positive, so that the critical velocity is increased by increasing 'k 2 . 
Thus in Ex. (1) k 2 was taken equal to |ft 2 , and the critical velocity 
obtained was J (774:0,) ; if, instead, we had taken ffi = a 2 , the critical 
velocity would have been iS J( 1 240ft).— Jan., 1904.] 

8. Character of the Fluctuations about Steady Motion, Mode in which the 

System Overturns. 

The character of the fluctuations about steady motion depends on 
the nature of the roots of the biquadratic (4ft), the expressions for the 
displacements 8u, §v, BO being evidently of the following forms :— 

For roots all real . ...................... cjgM + e 2 e^ + e 8 g V -f c±e K A 

For roots two real, two imaginary . . . Cjg'V + c^eM -f ye at cos (fit - c). 

For two pairs of imaginary roots ...... y^i* cos (/3ji ~~ e 2 ) 

+ y^eM cos (P4 - €2). 

The last form indicates two different sets of undulations of different 
lengths. Photographs of the paths of gliders taken by magnesium 
light distinctly show these two undulations, thus confirming our 
theory. 

An important further consequence is that a glider may perform 
undulations decreasing in amplitude, corresponding to a pair of complex 
roots of the biquadratic with their real part negative, but the motion 
may be unstable through the other roots having their real part positive, 
or one or both of them being real and positive. This indicates a real 
danger in experimenting with gliders. 

Stability may be broken either if a real root of the equation (4ft) 
changes from negative to positive, or if the real part of a pair of 
imaginary roots changes from negative to positive. 

The condition for the latter is that H = 0. whereas a real root 



1904.] The Longitudinal Stability of Aerial Gliders. 115 

changes sign if E = 0, If stability is broken by a fall of velocity the 

only quantities which can vanish are H and 0, of these II is more 

difficult to make positive than 0, and hence it appears that the most 

likely way for a glider to overturn in general is by commencing with 

a series of oscillations of increasing amplitude. This again agrees with 

experiments. 

9, Conclusions, 

1. For a glider or other body moving in a vertical plane in a 
resisting medium of any kind whatever, the small oscillations about a 
state of uniform rectilinear motion are determined by an equation of 
the fourth degree, so that the conditions for stable steady motion are 
those obtained by Routh. 

2. The coefficients in the period equation involve, in addition to the 
ordinary dynamical constants, nine quantities 'K u . . , G#, which, when 
referred to rectangular axes fixed In the body, represent the differential 
coefficients of the forces and couple due to the aerial resistances with 
respect to its translatory and rotatory velocity components. 

3. In the case of a system of aeroplanes these nine quantities can 
be expressed for the separate planes in terms of /'(a) and <£'(a), 
where f(a) and <£(a) are functions determining the resultant thrust, 
and the position of the centre of pressure when the direction of the 
relative motion of the air makes an angle a with the plane. These 
functions have been tabulated for certain different forms of surfaces, 
but further data are greatly needed. 

4. The longitudinal stability of the gliders is thus seen to be capable of 
mathematical investigation, and it is of paramount importance that the 
present methods should be practically applied to any aerial machines 
that may be designed or constructed before any actual glides are 
attempted. 

5. The methods of calculation are exemplified by numerical deter- 
minations of the criterion of stability in the cases of a single plane 
lamina and a pair of planes one behind the other. Most of the calcula- 
tions have been performed for an angle of gliding of 10° with the 
horizon, and it has been necessary to assume arbitrary values for the 
moment of inertia of the lamina. 

6. The condition that any steady linear motion may be stable in all 
these cases assumes the form V 2 > Jca, where a is a constant depending 
on the linear dimensions of the glider, and k is a constant depending 
on Its shape, the angle of gliding and the law of aerial resistance. 

7. For a pair of narrow slats, in which the variations in the positions 
of the centres of pressure of each are neglected, certain coefficients of 
stability vanish. If the planes are square so that the displacements 
of the centres of pressure are not neglected, the system is less stable 
than a single plane of breadth equal to one of the squares. 

8. By inclining the planes at a small angle to each other the stability 
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is much increased. On the other hand, if they are made to slope away 
from each other the glider becomes unstable. 

9. Two square planes of equal size placed one behind the other at a 
small angle are less stable in the examples considered than a square 
equal in size to one of the two, but more stable than a single square 
whose side is equal to the total length of the glider formed by 
the pair. 

10. A pair of unequal squares of which the smaller forms a rudder 
are more stable, in the examples considered, whe'n gliding with the 
rudder behind than with the rudder in front. 

11. In general, the stability is increased by making the moment of 
inertia as small as possible. 



" Cultural Experiments with ' Biologic Forms ' of the Erysiphacem." 
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(Abstract.) 

: In the introductory remarks the author points out that through 
specialisation of parasitism " biologic forms " have been evolved in the 
Erysiphacece which, both in their conidial (asexual) stage and ascigerous 
(sexual) stage, show specialised and restricted powers of infection. 
The powers of infection, characteristic of each " biologic form," are 
-under normal conditions sharply defined and fixed, and hitherto the 
result of the experiments of numerous investigators — both in regard 
to the present group of fungi and to the Uredinece, where the same 
specialisation of parasitism occurs — has been the accumulation of 
evidence tending to emphasise the immutability of " biologic forms." 

The second part of the paper gives the result of cultural experi- 
ments with " biologic forms n of Erysiphe Graminis DC, carried out 
during the past summer in the Cambridge University Botanical 
Laboratory. It has been found that under certain methods of culture, 
in which the vitality of the host-leaf is interfered with, the restricted powers of 
infection, characteristic of " biologic forms ," break down. 

In the first method of culture adopted, the leaf, which was either 
attached to a growing plant, or removed and placed in a damp 
chamber, was injured by the removal of a minute piece of leaf-tissue. 
In this operation the epidermal cells on one surface, and all or most 
of the mesophyll tissue, Avere removed at the cut place, but the 
epidermal cells on the other surface (opposite the cut) were left un- 
injured. Conidia were sown on the cuticular surface of the uninjured 
epidermal cells over the cut. In a few experiments the conidia were 



